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Abstract 

We illustrate the current status of heavy quark physics on the lattice. Special 
emphasis is paid to the question of systematic uncertainties and to the connec- 
tion of lattice computations to continuum physics. Latest results are presented 
and discussed with respect to the progress in methods, statistical accuracy and 
reliability. 



1 Introduction 

The role of lattice QCD in the field of heavy quark physics is twofold. 



First of all it provides experimental and phenomenological physicists with reliable 
predictions of those non-perturbative QCD parts, which are in nature almost inevitably 
connected to most of the important electroweak processes. For example, the Cabbibo- 
Kobayashi-Maskawa matrix elements for the weak decay of the b quark cannot be 
extracted from experimentally measured B meson decays without distinct knowledge 
about QCD bound state properties of heavy light systems, i.e. the form factors. 

Secondly, in the heavy quark regime lattice QCD can successfully describe "pure QCD" 
quantities like quarkonia splittings and the heavy quark potential, which, in turn, can 
be used to determine the strong running coupling constant a^. 

*Invited talk given at HEAVY QUARK PHYSICS AT FIXED TARGET, St. Goar, Oct. 3rd, 
1996. 
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In order to judge on the quality of results from lattice QCD it is necessary to under- 
stand roughly how the lattice method works in practice, what are its advantages and 
shortcomings, where systematic uncertainties originate from and how these uncertain- 
ties can be controlled. 

In principle, the lattice procedure consists in solving the QCD path integral on a 
(Euclidean) space time lattice. The connection with continuum physics is made in 
the limit, where (a) the lattice constant a goes to zero, and (b) the lattice volume V 
goes to infinity. If it were possible to achieve this limit unambiguously, QCD would 
be solved. 

In practice, however, the path integral can be evaluated only for finite values of a and 
V. This is done numerically with the help of Monte Carlo methods, which introduce a 
statistical uncertainty on the results. In order to achieve the continuum limit, extrap- 
olations in a and V are necessary. On top of this, most of the lattice calculations in the 
past have been performed in the so called quenched approximation, where (roughly 
speaking) internal fermion loops are neglected. 

Progress in lattice QCD therefore has to be measured in terms of 

(a) Statistical significance 

(b) Reliability of the extrapolation a — > 0, F — > oo 

(c) Ability to include internal fermion loops. 

Especially in the case of heavy-light systems, item (b) confronts lattice QCD with a 
serious problem. The inverse of the lattice constant can be viewed as a (gauge 
invariant) cutoff to all lattice observables. With todays computer facilities, cutoffs up 
to ~ 3.5GeV can be achieved for suitably large lattice volumes. If one would try 
to calculate B-meson properties , with Mb — 5GeV, directly on such a lattice, severe 
finite cutoff effects would prevent a reliable extrapolation a — > 0. 

An obvious way to alleviate this problem is of course to extrapolate lattice data in 
mass from far below the cutoff into the region of heavy masses, e.g. to Mb- However, 
as the functional dependence of the lattice data on the heavy quark mass is in general 
not exactly known, this procedure introduces additional systematic errors, which need 
to be controlled. 

In order to overcome this unsatisfactory situation, several methods have been devel- 
oped and improved over the recent years. These methods can be classified into three 
different groups: a priori methods, phenomenological methods, and effective methods. 
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A priori methods trie to improve on the discretization of the QCD action in such a way, 
that cutoff effects are reduced. Those "improved actions" converge faster to the correct 
continuum form, leaving the QCD physics unchanged. The most prominent example[] 
in this group is the Sheikholeslami-Wohlert action [Q, which is often called clover ac- 
tion. In contrast to the standard Wilson action^], finite cutoff effects proportional 
to a are absent, at least on the classical level. It turns out, however, that quantum 
fluctuations can re- introduce 0(a) cutoff effects. This problem is currently tackled by 
a proper adjustment of the coefficient of the (additional) clover term. Lattice actions 
which include this adjustment are traded under the names "tadpole improved clover 
action"!^ or "non-perturbatively improved clover action" [§. 

In principle it should be possible to formulate a "perfect" lattice action, which is free of 
any cutoff effect. The construction of such an action is currently under consideration 0. 

The basic idea of phenomenological methods is, to reduce finite cutoff effects by proper 
adjustment of lattice observables. This is achieved by a change in the normalization 
of the lattice quark propagator (LMK I) and by a redefinition of the particle masses 
(LMK 11)0. As the modifications are not implemented on a fundamental level, it is 
not clear whether this method leads to a general improvement of all lattice observables. 

Effective methods have been designed for heavy quarks on the lattice. The idea is 
to remove the largest scale, i.e. the heavy quark mass, from the lattice action. If 
all remaining scales are small compared to the lattice cutoff, finite a effects should 
be reduced substantially. Lattice implementations of effective methods have been 
developed for (a) the static approximation |T^, i.e. the zeroth order of heavy quark 
effective theory, and (b) the non-relativistic QCD[jll| (NRQCD). The latter approach 
starts from a non-relativistic approximation of QCD and includes, similar to the well 
known Fouldy-Wouthusen transformation, successively relativistic corrections in form 
of a 1/Mq expansion. Clearly, the range of validity of such methods is limited to the 
region of (very) heavy quarks. 

In view of the variety of all these methods, whose merits in some cases have not been 
fully proven yet, one could argue that lattice QCD has lost its status of being an ab 
initio method. In the following we will try to convince you of the opposite. 

The backbone of large scale lattice calculations is always a powerful computer. The 
strong increase in sustained computer speed over the recent years allowed for a sub- 
stantial improvement on statistical significance of lattice data as well as for a variation 
of the lattice constant and the lattice volume over a much larger range. With such 
a tool in hands, one can indeed check on the benefit of the different methods and 

^We comment here only on the fermionic part of the QCD action. For improvements on the gluonic 
part see refs.[^ ||, ||. 



3 



calculate a reliable estimate of the size of systematic uncertainties. For example, one 
can measure the size of finite cutoff effects by performing a series of lattice simulations 
with different lattice constants a. In that view, the computer can be compared with 
the accelerator and the detector of a large experiment: the quality of the device has 
an decisive influence on the quality of the results. 

In the following the actual status of the "lattice experiments" within the field of 
heavy quark physics is reviewed in form of selected examples. We will discuss the 
decay constant of the B meson, fs-, the semileptonic decays of B mesons, and the 
determination of from quarkonium splittings and the heavy quark potential. 



2 Status of Jb 



2.1 Strategy 



The decay constant of the B meson, is defined in the space time continuum by 

(0|Ao|i?)^°"* = mB/B , (1) 

where Aq denotes the zeroth component of the axial current and is the B meson 
mass. Thus, can be determined once the matrix element on the left hand side of 
eq.(|l|) is known. 

On the lattice, one calculates (0|y4o|-B)''^**, which is related to its pendant in the con- 
tinuum by 

(0|Ao|B)"°"* = Z^(0|Ao|5)'"" . (2) 

The renormalization constant Za accounts for the non- conservation of the axial cur- 
rent. Both, (0|y4o|-B)'"** and depend on the lattice spacing a, and much work has 
been devoted to the problem of choosing Za such that it just cancels the finite cutoff 



effects of the above product [12, 13 1. As a satisfactory solution to this problem has not 



been found yet, an extrapolation a ^ is still necessary. 



Due to limitations in computer time - and computer memory size - the B meson 
mass is beyond currently attainable lattice cutoffs. Therefore, one calculates fps with 
meson masses below the cutoff and in the static limit. An 1/mps expansion then 
interpolates between these results, yielding at a given value of the lattice cutoff. 
The quality of such an interpolation is demonstrated in fig.0. 

This procedure has to be repeated for several cutoff values, and finally the continuum 
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Figure 1: MILC collaboration^^ : fqqyMqq vs. I/Mq^ at a finite value of the cutoff 
(a~l ^ 3.2GeV). The solid line is a quadratic fit to the octagons ( "heavier heavies" 
+ static); the dashed line is a quadratic fit to the crosses ("lighter heavies" +static). 
The difference of the two curves in the heavy quark region illustrates the quality of the 
interpolation. 



extrapolation a — has to be performed^. 

As an alternative to this standard procedure, one can in principle use the NRQCD 
effective method. Unfortunately, NRQCD is a non-renormalizable theory and one 
cannot send a ^ 0. Therefore, one needs good control over discretization errors of the 
NRQCD action. In the end it is necessary to verify, that the results from the standard 
method and from NRQCD are consistent. 



2.2 Results 

Three years ago, the first calculation, which includes the a — > extrapolation^, has 



been performed by the PSI-WUPP collaboration! 1 In this simulation, the cutoff and 



the lattice volume were varied in the ranges l.SGeV < a ^ < 3.0GeV and l.Ofm < 



continuum extrapolation includes in general also the limit V ^ oo. As the finite volume effects 
are not crucial for heavy quarks, we will not discuss this point in detail. 
■^Results at fixed cutoff have been published by |l^, [l^ . 
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V < 1.5fm respectively. The continuum extrapolation yielded Jb = 180(50)MeV. 

This year, the lattice determination of Jb has been improved by large scale simulations 
of the JLQCD group[l^ and the MILC collaboration ||19||. Preliminary results of a 
NRQCD calculation of fs have been published by the SGO collaboration . 
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Figure 2: JLQCD collaboration!^^ : Jb^ as a function of the lattice spacing a. Dif- 
ferent symbols refer to different choices of Z^: naive (open circles), KLM (open dia- 
monds), KLM improved (crosses), mixed KLM (filled circles). 

JLQCD has pushed forward in the size of cutoff and volume, namely l.SGeV < < 
3.4GeV and 1.7fm <V< 2.0fm. An important ingredient of their analysis is a careful 
study of the influence of different choices of Za on the cutoff dependence of /b- The 
result is displayed in fig.^. It turns out that, although none of the choices of Za 
removes the cutoff dependence completely, the continuum values agree within errors. 
The remaining systematic uncertainty due to Za can be read off the spread of results 
at a = 0. The final JLQCD results read 



fB = 179(ll)12^MeV /b. 



197(7)1^5MeV 



The first error accounts for statistical, the second for systematic uncertainties. It 
is very encouraging to see the PSI-WUPP result being consolidated by the - more 
accurate - JLQCD calculation. 



The MILC collaboration has extended the cutoff even further. Their largest value is 
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Figure 3: MILC collaboration /[T^j : fs^ cls a function of the lattice spacing a. Octagons 
refer to quenched data; crosses to full QCD, nj = 2. The solid line is a fit to all 
quenched points (conf. level = 0.66); the dashed line is a constant fit to the three 
quenched points with a < Q.^GeV (conf. level = 0.76). The extrapolated values at 
a = are indicated as bursts. 



~ 4.5GeV. Fig.^] shows the final data as well as the extrapolation to a = 0. For 
the first time in lattice QCD, MILC has been able to include full QCD simulations - 
with rif = 2 dynamical fiavors - into the determination of fs- The corresponding data 
is represented by crosses in fig.|^. For small a the full QCD data seem to enhance the 
value of fs- However, the statistical uncertainty is too large to draw a firm conclusion. 
Instead, MILC quotes the quenched result and includes the effect of dynamical fermions 
ClS db systematic uncertainty 

fe = 166(ll)(28)(14)MeV fe, = 181(10)(36)(18)Me\/ . 

The first error includes statistical errors and systematic effects of changing fitting 
ranges, the second other systematic errors within the quenched approximation, and the 
third accounts for quenching effects. Within errors, the MILC result is well compatible 
with PSI-WUPP and JLQCD. 
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The SGO collaboration has analyzed Uf = 2 full QCD configurations, using NRQCD 
for the heavy quark. They quote a (preliminary) value of fs — ISOMeV. 

In summary, the "old" value fs = 180(50)MeV has been consolidated by more ad- 
vanced lattice simulations. However, much work remains to be done in order to reduce 
statistical and systematic uncertainties. A major goal of the next years will be the 
inclusion of dynamical fermion loops. 



3 Semileptonic Decays of B Mesons 

Semileptonic decay amplitudes can be described as a product of a (perturbatively 
accessible) weak interaction term and perturbatively not accessible QCD part. The 
latter is commonly parameterized by a set of form factors 

(PSVmIPS) = Expr{/+(g2),r(g2)} (3) 
(V|J,|PS) = Expr{v(g2),A(g2),A2(g^),A(g2)} . 

PS, PS' and V are momentum dependent pseudo-scalar and vector meson states, and 
is the (axial-) vector current. Clearly, the contribution from lattice QCD is a deter- 
mination of the form factors by a calculation of the left hand side of eq. (^ . Compared 
to the calculation of this is not an easy task, as both, the mass dependence and 
the behavior has to be determined. Due to limited statistical accuracy, it is not 
yet possible to extract the functional behaviors unambiguously from the data. Thus, 
systematic uncertainties due to the choice of inter- and extrapolation methods have 
to be estimated and finally included into the results. 



3.1 B ^ TT, p 

Considerable progress with respect to the reliability of inter- and extrapolations has 
been achieved within the last two years. "Old" methods[^, 23|, simply relied on the 



validity of the pole dominance hypothesis and on HQET. As the latter is justified only 
in the region ~ g^ax had to operate at very large values of g^. 

As current lattice investigations substantially improved on the statistical quality of 
data, it has become possible to work at more moderate values and to check on the 
infiuence of the various assumptions. 



The Wuppertal group [p4[ has performed its analysis at small values of g and for meson 



masses up to Md. In this region the form factors can be determined quite reliably. 
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The extrapolation in mass was done considering various parameterizations for the mass 
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Wuppertal group^^J: Mass dependence of form factors F{q^ = 0) for the 
> vr and B ^ p. Crosses denote the data. The solid line refers to a fit to 



(a), the dashed line to F = -^7^^ (02 + Af^s) dotted line 

('^3 ~^ Ivies') '^^^ extrapolated results are depicted at l/mB[GeV~^] 



Figure 4: 

decays B 
F = ai^ 

to F = a/Mp^s 

as filled circles (method a), upper triangles (method b), and stars (method c). 

dependence, as shown in fig.||. Unfortunately the data is still not precise enough to 
discriminate between the Ansatze. However, the systematic uncertainty due to the 
extrapolation in mass can be estimated from the variation of results with respect to 
the various Ansatze. 



An analysis which is more guided by the Heavy Quark Effective Theory (HQET) has 
been performed by the UKQCDpSf collaboration. In order to extract the maximal 
information from the data, they work at both, small and large values of g^, combining 
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Table 1: Form factors at q"^ = for the semileptonic decays B vr(/z/), B p(^z^) 



/+(0) 



/o(0) 



y(o) 



AM 



^2(0) 



0.24+^ 
0.50(14)t^ 
0.29(6) 
0.35(8) 
0.28(14) 
0.33(17) 



3 



0.27+^ ± 3 
0.16(4)1?2 

0.29(16) 
0.24(12) 
0.24(6) 
0.21(5) 



+y +4 

-6 -5 



UKQCD 
Wupp[|| 

APE„ra 



APEfop 
ELa,[|22| 



ELCfep2| 



0.24 
0.20(3) 



0.61(23)1^ 
0.45(22) 
0.53(31) 
0.37(14) 
0.40(16) 



0.28 



0.72(35)1^° 
0.24(56) 
0.27(80) 
0.39(24) 
0.47(28) 



Table 2: Branching ratios (Vo) for exclusive semileptonic B decays with a D or D* 
final state 





B^ D+l-y 


B D*+l-u 


UKQCD 
WUPP 
CLEO I 
CLEO II 
ARGUS 
ALEPH 


1.5lt±0.3 
2.1(1.4) 
1.8 ±0.6 ±0.3 

2.1 ±0.7 ±0.6 


4.8t^±0.5 

4.1 ±0.5 ±0.7 
4.49 ±0.32 ±0.39 

4.7 ±0.6 ±0.6 
5.36 ±0.50 ±0.76 



the results after extrapolation to Mb- Although this method assumes the validity of 
HQET at moderate meson masses, e.g. at the D meson mass, it leads to quite definite 
conclusions on the dependence of form factors at the B meson mass. 

Table |1| shows the results of the different groups^. It turns out that the uncertainties 
are still too large to try a reliable extrapolation a — > 0. 



3.2 B^D,D* 

In order to estimate the systematic effects due to the various inter- and extrapolation 
methods, the Wuppertal group has applied the same type of analysis as for the decays 
B ^7r,p. 

In contrast, the UKQCD collaboration has performed its analysis completely within 
the framework of HQET. One major conclusion of their work is, that non-perturbative 
"power corrections" to the form factors are small in the mass region of D and B mesons. 



In table |^ the results for the branching ratios are compared to the measurements ||27 
of ARGUS, CLEO and ALEPH. 



*UKQCD uses /+(0) = /o(0) as a constraint. 
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4 Quarkonia splittings and ag 



Pure heavy quark systems like quarkonia provide an ideal laboratory to study (gluonic) 
inter quark forces. Lattice QCD is well prepared to work in this lab with tools like 
NRQCD and improved actions. 



In order to demonstrate their quality, we compare in fig.^ the experimentally measured 
T spectrum with the lattice results [^]. Both, the NRQCD data and the results 
of the calculation with tadpole improved clover action are in good agreement with 
experiment, once dynamical fermions are included. This sets the stage to extract the 
strong coupling from the inter quark forces in the heavy quark regime. 
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Figure 5: T spectrum (left) and spin splittings. The data has been compiled by 
J. Shiqemitsuf^ . 



In principle, any lattice quantity which can be reliably expanded in powers of a can be 
used to define the strong coupling on the lattice. In quarkonia studies one commonlyQ 
works with ap[^, which is defined by 



47r 



3.4 



a 



InWi^i = —ap — [1 - (1.19 + 0.07n/)ap] 



(4) 



■^Another ingenious choice is agfp9[|, which is used in connection with the Schrodinger functional 
technique. Unfortunately, full QCD results for asp are not available yet. 
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Wi^i is the 1x1 Wilson loop, which has to be measured on the lattice. In order 
to determine the momentum scale Q = ^ of ap, the lattice cutoff needs to be 
known. This is the point, where quarkonium physics enters: can be extracted by 
comparing the measured IS — 2S or S — P splittings with the lattice data. Finally, 
extrapolations in the number rif of dynamical flavours and in the mass of dynamical 
quarks rriq have to be carried out, n/ = 0, 2 — > ri/ = 3, nig TTT^eff = {rnu+md+ms)/3. 
As a result of this procedure, the NRQCD collaboration quotes pif] 

r 0.1948(29) (11) (37) S - P splitting 

«?'='^[8.2GeV] = 

[ 0.1962(41) (08) (40) IS - 2S splitting 

where the first error accounts for statistical uncertainties, the second for discretization 
effects and the third for uncertainties due to the extrapolations. To compare with 
other determinations of the strong coupling it is advantageous to convert ap to aj^. 



Unfortunately, the corresponding two loop conversion formula has been calculated |32 
hitherto only for n/ = 0. This leads to an additional systematic uncertainty, which 
has to be taken into account. At the Mzo mass, the NRQCD collaboration finds p8[] 

r 0.1175(11)(13)(19) S - P splitting 

[ 0.1180(14)(14)(19) IS - 2S splitting 

Note that the last error, which accounts for the uncertainty in the conversion, domi- 
nates. 



5 Summary and Conclusions 

We have discussed the importance of recent developments and results in heavy quark 
physics on the lattice. Considerable progress has been achieved with respect to the 
reliability of the continuum extrapolation and to the inclusion of internal fermion 
loops. This is mainly due to the dramatic improvement in computer performance. 

In order to reduce discretization effects, a variety of promising new ideas and methods 
has been developed. With the advent of fast parallel computers, the benefit of these 
methods can be tested reliably in numerical simulations. Once this has been done, 
they will be used for precise calculations at or even beyond the b quark mass. 
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